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Lepton flavor violating Higgs boson decays are studied within the context of seesaw models 
with Majorana massive neutrinos. Two models are considered: The SM-seesaw, with the Standard 
Model Particle content plus three right handed neutrinos, and the MSSM-seesaw, with the Minimal 
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^ ' Supersymmetric Standard Model particle content plus three right handed neutrinos and their su- 

(N ■ 

■ persymmetric partners. The widths for these decays are derived from a full one-loop diagrammatic 

m 

computation in both models, and they are analyzed numerically in terms of the seesaw parameters, 



namely, the Dirac and Majorana mass matrices. Several possible scenarios for these mass matrices 



' that are compatible with neutrino data are considered. In the SM-seesaw case, very small branch- 

^ • ing ratios are found for all studied scenarios. These ratios are explained as a consequence of the 

r~| , decoupling behaviour of the heavy right handed neutrinos. In contrast, in the MSSM-seesaw case. 



sizable branching ratios are found for some of the leptonic flavor violating decays of the MSSM 
neutral Higgs bosons and for some choices of the seesaw matrices and MSSM parameters. The 
relevance of the two competing sources of lepton flavor changing interactions in the MSSM-seesaw 
case is also discussed. The non-decoupling behaviour of the supersymmetric particles contributing 
in the loop-diagrams is finally shown. 
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I. INTRODUCTION 

The present strong evidence for lepton flavor changing neutrino oscillations W in solar 
and atmospheric neutrino data, as well as in the KamLAND reactor experiment, implies the 
existence of non-zero masses for the light neutrinos, and provides the flrst experimental clue 
for physics beyond the Standard Model (SM). The experimentally suggested smallness of the 
neutrino masses can be explained in a very simple and elegant way by the seesaw mechanism 
of neutrino mass generation |2l]. This mechanism requires the introduction of heavy right- 
handed (RH) Majorana neutrinos which are singlet under the SM gauge symmetry group and 
whose Majorana masses, m^v/^ , can therefore be much higher than the SM particle masses. In 
this context, the smallness of the light left-handed (LH) neutrino masses appears naturally 
due to the induced large suppression by the ratio of the two very distant mass scales that are 
involved in the seesaw mass matrices, the Majorana matrix tum and the Dirac matrix 171^. 
The latter is generated after electroweak symmetry breaking by mn = < H >, where Y^, 
is the Yukawa matrix for couplings between the RH and LH neutrinos, and < H >= v = 174 
GeV is the SM Higgs boson vacuum expectation value. For the one generation case, and 
assuming a Yukawa coupling of order one, the suggested small neutrino mass value signals 
towards a new phyics mass scale of the order of tum ~ 10^^ GeV, but the pattern and size 
of the seesaw mass parameters can vary much respect to this in the most general case of 
three generations. 

Another appealing feature of the seesaw models is that the RH Majorana neutrinos can 
succesfully generate, through their CP-violating decays and via leptogenesis, the observed 
baryon asymmetry of the Universe. There is, however, one negative aspect in the standard 
version of the seesaw models. It is that the presence of the two distant mass scales can lead to 
a severe hierarchy problem, and this requires the introduction of supersymmetry (SUSY) to 
be solved. In the SUSY-seesaw models, the hierarchy between ttlm and the electroweak scale 
is stabilized by the new contributions of the SUSY partners of the RH and LH neutrinos. 
Thus, the SUSY-seesaw models, and particularly the simplest version given by the Minimal 
Supersymmetric Standard Model (MSSM), are becoming more popular. 

One of the most interesting features of the SUSY-seesaw models is the associated rich 
phenomenology due to the occurrence of lepton flavor violating (LFV) processes. Whereas 
in the standard (non-SUSY) seesaw models the ratios of LFV processes are small due to 
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the smallness of the hght neutrino masses, in the SUSY-seesaw models these can be large 
due to an important additional source of lepton flavor mixing in the soft-SUSY-breaking 
terms. Even in the scenarios with universal soft-SUSY-breaking parameters at the large 
energy scale associated to the SUSY breaking Mx (which could be the Planck mass, the 
SUSY-GUT mass or something else, but always well above m^), the running from this 
scale down to tum induces, via the neutrino Yukawa couplings Y^, large lepton flavor mixing 
in the slepton soft masses, and provides the so-called slepton-lepton misalignment, which 
in turn generates non-diagonal lepton flavor interactions. These interactions can induce 
sizable ratios in several LFV processes with SM charged leptons in the external legs, which 
are actually being tested experimentally with high precision and therefore provide a very 
interesting window to look for indirect SUSY signals. 

In addition to the previously mentioned radiatively induced LFV effect, there is another 
source of radiative LFV effects in the MSSM-seesaw, which is generated from the neutrino 
sector mixing. Namely, the off-diagonal entries of the Maki-Nakagawa-Sakata (MNS) matrix, 
Umns 0, produce lepton flavor non-diagonal interactions involving neutrinos which could 
also generate, via loops, important contributions to the LFV processes with SM leptons 
in the external legs. In the MSSM-seesaw, these LFV effects from neutrino mixing could 
indeed compete in some processes with the other source from slepton-lepton misalignment. 
This competitivity of the two flavor changing sources has indeed been noticed previously 
in the squark sector, concretely in the SUSY electroweak one-loop contributions to flavor 
changing Higgs boson decays, bs, sb 0. Interestingly, the lepton sector could manifest 

more evidently these competing effects than the quark sector, since the MNS matrix being 
prefered by neutrino data is clearly non-diagonal, whereas the CKM matrix is close to the 
identity matrix. 

Among the various LFV processes that have been considered in the literature, the most 
fruitfull ones are the radiative — 67, r — /i7,_a.iid r cy decays 

mm 

, since their 



These usually provide the most 



restrictive experimental bounds on the MSSM-seesaw parameters. Another interesting LFV 
decays include r rare decays, Z boson decays Z 1^1^ and Higgs boson decays H Ikim, 
with k ^ m. 

We are interested here in the LFV Higgs boson Decays (LFVHD), H — > Tfi,Te,fie, 
and the branching ratios that can be generated in the context of the seesaw models with 
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parameters being compatible with the neutrino data and the most relevant data of r and 
fi radiative decays. In this work we consider both versions of the seesaw mechanism for 
neutrino mass generation: the SM-seesaw and the MSSM-seesaw. In the first one, the SM 
particle content is enlarged by three singlet RH Majorana neutrinos, and in the second one, 
the MSSM particle content is enlarged by three singlet RH Majorana neutrinos and their 
corresponding SUSY partners. We present here a complete one-loop computation of the 
LFVHD widths in both seesaw scenarios and analyze the size of the associated branching 
ratios in terms of the seesaw parameters. In order to make contact between the seesaw mass 
matrices and the experimental neutrino data, we use the general parametrization introduced 
in ref. {gJ, where the Yukawa neutrino couplings, and therefore m^), are expressed in terms 
of the three physical light neutrino masses, m,^., the three physical heavy neutrino masses, 
rriN^, the Umns matrix, and a general complex 3x3 orthogonal matrix R. 



m. There, it 



The LFVHD within the SM-seesaw was studied some time ago in ref. 
was considered a very specific seesaw scenario where all the light neutrinos were exactly 
massless at tree level [l^ and the Dirac mass was taken very large. The conclusion was 
that branching ratios for Hsm Tfi decays as large as 10~^ — 10^^ can be achived for 
^HsM — GeV. Besides, these ratios were found to grow with the heavy neutrino masses, 
and this growing suggested a non-decoupling behaviour of the heavy neutrinos. The crutial 
assumption for that behaviour was to take very large values, which implied very strong 
neutrino Yukawa couplings. Here, we have recalculated the LFVHD rates in the complete 
one-loop diagrammatic approach and updated the numerical estimates in the light of the 
recent neutrino data, by fixing the input parameters mj,^ and Umns to their data prefered 
values. We then make our estimates for specific choices of the input unknown parameters, 
rriNi ^^"^ ^^"^ P^y special atention to those which generate succesfull baryon asymmetry. 
For all the studied cases we find ratios that are many orders of magnitude smaller than in 
ref. The reason for these small ratios is because the heavy RH neutrinos do indeed 

decouple in our SM-seesaw scenario and our assumptions for the seesaw parameters do not 
imply large rriD values. For completeness, we also include in this part, an estimate of the 
LFVHD rates in the case of Dirac massive neutrinos and compare them with the previous 
estimate of ref. jl4|. 



The second part of this work concerns with the evaluation of the LFVHD ratios in the 
context of the MSSM-seesaw. Concretely, h^,H^,A^ — > Iklmi with k ^ m. The subject 
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of LFVHD being generated from loops of SUSY particles has been considered previously 
in refs. and 3|- In ref. Q] it was analyzed a specific SUSY-SU(5) scenario where 



rom the running of the trilinear 
was not in the context of the 



the slepton-lepton misalignment was generated exclusive 
A-terms. On the other hand, the computation of ref. 
MSSM-seesaw but in a more generic scenario for slepton-lepton misalignement. Besides, 
in ref. Q|, the effective lagrangian approach that is valid for large tan/3 values and large 
SUSY mass values is used. We present here instead, a complete one-loop computation 
in the MSSM-seesaw context and do not rely on any of the above approximations. We 
include in the computation both sources of lepton flavor violating interactions, the slepton- 
lepton misalignment and the neutrino mixing via Umns- The slepton-lepton misalignment 
is generated, as it is usual in the seesaw models, by the renormalization group running of the 
slepton soft parameters from the high energies Mx down to mAf. The diagonalization of the 
generated slepton mass matrix is then performed. That is, we do not rely on either the mass 
insertion approximation or the large tan/3 effective lagrangian approach and, therefore, our 
results are valid for all tan (3 values and all soft-SUSY-breaking mass values. We will explore 
the size of the branching ratios for the Higgs decays as a function of the relevant MSSM 
parameters, which within the context of mSUGRA are the universal soft masses Mq, M1/2 
and tan/3, and of the relevant seesaw parameters, which are and the R matrix. We will 
analyze in parallel the branching ratios for the Ij —>■ k'j decays as a function of the same 
parameters. The requirement of compatibility with the present data on Ij —>■ lij decays, 
mainly fi ^ e'y and t ^ fi'y, will provide us with the maximum allowed ratios for the Higgs 
decays. We will also study the behaviour of the LFVHD widths in the limit of very heavy 
SUSY masses and will find that the sleptons, sneutrinos, charginos and neutralinos do not 
decouple in this observable. For large SUSY masses, large tan (3 and particular choices of 
the see-saw parameters we will find agreement with the numerical results of ref. jisl ]. 

The paper is organized as follows. Section II contains a short summary of the mass 
parameters and mixings in the neutrino sector of the seesaw models. The relation between 
these parameters and the experimental neutrino masses and mixings is also included. Section 
III is devoted to the computation and analysis of the LFVHD rates in the context of the 
SM-seesaw. The decoupling behaviour of the heavy neutrinos is also studied in this section. 
Section IV starts by presenting the two sources of LFV interactions in the context of the 
MSSM-seesaw. Next the computation and analysis of the LFVHD rates within that context 
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is included. The non-decoupling behaviour of the SUSY particles in the LFVHD widths is 
studied at the end of this section. Section V is devoted to the conclusions. 



II. NEUTRINO MASSES AND MIXINGS IN THE SEESAW MODELS 



In this section and in order to fix our notation we briefly review the mass parameters and 
mixings in the neutrino sector of the seesaw models and relate them to the physical light 
neutrino masses and neutrino mixing angles which are extracted from neutrino data. We 
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3 



for SM-seesaw and 



follow closely here and in the next sections the notation of refs. 
ref. joj for SUSY-seesaw and for the connection with neutrino data. 

We start with the Yukawa-sector of the SM-seesaw that contains the three LH SM neu- 
trinos i/^j and three extra RH massive neutrinos z/^ ^, whose Yukawa interactions provide, 
after spontaneous electroweak symmetry breaking, together with the right handed neutrino 
masses, the following mass Lagrangian containing the Dirac and Majorana mass terms, 

-^:;a.. = ^K,Km^|^^"'^; I +h.c., (i) 

where, 

I • (2) 

mj^ vfiM 

Here m£, is the 3x3 Dirac mass matrix that is related to the 3x3 Yukawa coupling matrix 
and the SM Higgs vacum expectation value, < H >= v = 174 GeV, by itld = < H >; 

and rriM is the 3x3 Majorana mass matrix for the RH massive neutrinos that is real, non 

singular and symmetric. 

The mass matrix M'^ is a 6 x 6 complex symmetric matrix that can be diagonalized by 

a 6 X 6 unitary matrix in the following way: 



U^T j^yjjy _ £j-v _ (^2a(7(m^j, m;^2, m,^3, m^r-^, m^Va, ?7^7V3)- (3) 

This gives 3 light Majorana neutrino mass eigenstates z/j, with masses m^j^ (i=l,2,3), and 
three heavy ones A^j, with masses mj^- (i=l,2,3), which are related to the weak eigenstates 
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via, 





j^.,, , u- 1 I and I I = :r I . (4) 

The seesaw mechanism for neutrino mass generation assumes a large separation between the 
two mass scales involved in and tum matrices. More specifically, we shall assume here 
that all matrix elements of mo are much smaller than those of itlm-, ttlo << mM, and we 
will perform an analytical expansion of all relevant interaction parameters and observables 
in power series of a matrix defined as, 

^ = niDml}. (5) 

In particular, the previous diagonalization of the mass matrix M'^ can be solved in power 
series of ^ jl6|. For simplicity, we choose to work here and in the rest of this paper, in a 
fiavor basis where the RH Majorana mass matrix, mjv/, and the charged lepton mass matrix, 
M', are fiavor diagonal. This means that all fiavor mixing of the light sector is included in 
the mixing matrix Umns- By working to the lowest orders of these power series expansions 
one finds, in the fiavor basis, the following neutrino 3x3 matrices. 



mj, = -moi^ + O^moi^) ^ -mDmj]^m]j (6) 
ruN = mM + 0{mDi) — mM- 

Here, is already diagonal, but is not yet diagonal. The rotation from this fiavor 
basis to the mass eigenstate basis is finally given by the MNS unitary matrix, Umns- Thus, 

m^t"''^ = UM^sm^UMNS = diag{m^^,m^^,m^^), (7) 
m^"^ = niN = diag{mNi,m,N^,mN^), 

and the diagonalization of M'^ in eqs. Q and (jHl) can be performed by the following unitary 
6x6 matrix: 
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As for the Umns matrix, we use the standard parametrization given by, 



U 



MNS 



— ■S12C23 — Ci2S23Si3e*'' C12C23 — Si2S23Si3e"' S23C13 
y S12S23 — 0120235136**^ — C12S23 — 5120235136**^ C23C13 J 



dm^(l,6*°,6*^). (9) 



where Cij = cos 9ij and Sij = sin Oij . 

Finally, in order to make contact with the experimental data, we use the method proposed 
in ref 0. It provides a simple way to reconstruct the Dirac mass matrix by starting with 
the physical light and heavy neutrino masses, the Umns niatrix, and a general complex and 
orthogonal matrix R. With our signs and matrix conventions this relation can be written 
as , 

ml = ^m';^''^'/'Rmt-^'/'UUs (10) 
where, R^R = 1. Thus, instead of proposing directly possible textures for m^, one pro- 
poses possible values for rriN-^ ,mN^ ,171^3 and R, and sets m^-^ , m^^ , m^g and Umns to their 
suggested values from the experimental data. Notice that any hypothesis for R different 
from the unit matrix will lead to an additional lepton flavor mixing between the LH and RH 
neutrino sectors. Notice also that the previous relation holds at the energy scale mM, and to 
use it properly one must use the RGE to run the input experimental data m^**^^ and Umns 
from the low energies mw up to rriM- We have computed here these running effects by solv- 



ing the RGE in the one loop approximation 



17| . and we have considered the corresponding 



radiatively corrected neutrino masses and Umns matrix elements in our computation of the 
LFV rates. 

For the numerical estimates in this paper we will consider the following two plausible 
scenarios, at the low energies, being compatible with data: 

• Scenario A: quasi-degenerate light and degenerate heavy neutrinos, 

m^, = 0.2 eV , m^^ = m^, + , m^^ = m^, + , (11) 

rriNi = '^N2 = ^N-i = ^N 

• Scenario B: hierarchical light and hierarchical heavy neutrinos. 



6V" , = ^ Am2^, , m^.^ = v^Am^^^, (12) 
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In the two above scenarios, we will fix the input low energy data to the following values, 



ol 



0.008 eV, ^/A^ 



atm 



O.OSeV, di2 = 9soi = 30", 623 = 9, 



atm 



45°, = 0° and 



.^„.,.0(Seefo..sta,.e,.ef.Q). 

Regarding the matrix R, we will consider correspondingly one of the following three cases: 

• Case 

R = Rq = 1 (13) 
This is a reference case which is chosen here just because it is the simplest possibility. 



Case 1 



^ C2C3 — C1S3 — S1S2C3 S1S3 — C1S2C3 ^ 



R = Ri 



\ 



(14) 



C2S3 C1C3 - S1S2S3 -S1C3 - C1S2S3 

S2 S1C2 C1C2 

where Cj = cos6i, Si = sin^j and di, 62 and 63 are arbitrary complex angles. This 
parametrization was proposed in ref. Id] for the study of yU — > 67 decays. It has also 
been considered in ref. jll with specific values for the 6i angles to study the implications 
for baryogenesis in the case of hierarchical neutrinos. 



Case 2 



R = Ro = 



(15) 



A 



\ 



J 



(16) 



with = 1, and 

^ a b 

—a c 
-b -c 

Here, a, b, and c are three real parameters that are constrained by perturbativity of 
the Yukawa couplings. In particular, for a = 6 = c = and riii,^ ^ ^ ^ 0.2eV, it leads 
to k < (1.4,0.9,0.3) for niN,^^-, ^ (10^°, 10^^ 10^^) GeV respectively. This choice has 
been proven in ref. to provide successful baryon asymmetry via leptogenesis for the 
case of quasi-degenerate neutrinos. 

All that has been sumarized in this section applies to the MSSM-seesaw model as well. The 
only difference is that in this case, rriD is related to one of the two Higgs vacum expectation 
values hj mn = Yj^ < H2 > , where < H2 >= V2 = v sin /3. 
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III. LEPTON FLAVOUR VIOLATING HIGGS DECAYS IN THE SM-SEESAW 



In this section we compute and analyze the LFVHD widths, r{H — > Iklm) with k ^ 
m, within the context of the SM-seesaw with three RH neutrinos which has been shortly 
reviewed in the prev ious section. The branching ratios for these decays were studied some 



time ago in ref. jlj] in a particular scenario of the SM-seesaw where all the light neutrino 
masses were exactly zero at the tree level and the Dirac mass was taken very large. The small 
masses for the light neutrinos were then generated by the one-loop electroweak radiative 
corrections ^ . Besides, the numerical results of ref. for the case of 3 generations were 
obtained under particular approximations for the evaluation of the one-loop diagrams, as 
for instance, m;^ = 0, << 1, -^^^ << I5 the latter being nowadays clearly not a very 
good approximation given the present lower bound on mn from LEP of about 115 GeV. 
Here we have recomputed these LFVHD widths for all the channels and we have included 
all contributing one-loop diagrams without assuming any approximation. One of the main 
points of our analysis is the update of the numerical results by taking into account the present 
experimental neutrino data. In addition we will reanalyze the behaviour of these partial 
widths with the large heavy neutrino masses. In contrast to ref. where the important 
enhancement of the widths found with m^r- suggested a non-decoupling behaviour of the 
heavy neutrinos, we will find instead a clear decoupling behaviour. The main difference 
between ref. [l2[ and us is the assumption on m^), which in their work was taken very 
large (therefore, leading to strong neutrino Yukawa couplings), whereas in our case the 
assumptions for the input 171^^ and R, do not lead to large m^. 

We start by writting down the interactions that are relevant for the computation of 
the LFVHD widths to one-loop and in the mass eigenstate basis. Denoting the Majorana 
neutrino mass eigenstates collectively by (i.e., rii = z/j for i=l,2,3 and rii = A^j_3 for 
i=4,5,6), the relevant interactions between and W^, H and the Goldstone bosons can 
be written respectively as follows. 



^int = -r-^Hrii [(m„^ + m„.)i?e(C„^„ ) + i-f^irun, - m„.)/m(C„.„ )] rij, 

Ltt = ^^G^~ [hBkn, {m^PL - mn,PR)n,] + h.c (17) 
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where, the couphng factors Bi^n^ (i=l,2,3, j=l,..6) and Cmuj (i,j=lv6) are defined in 
terms of the If^ matrix of eq. 01 by: 



Bl.n, = U^* (18) 

3 

Cn.n, = ^l^^l] (19) 



k=\ 



Notice that our particular choice of diagonal charged leptons in the fiavor basis is equiv- 



alent to assume a V matrix in [1^ equal to the identity matrix. The expansions of these 
coupling matrices in power series of the ^ matrix can be easily derived from the expansion of 
in eq. |H1 For brevity, we omitt here the indices and use a self-explanatory short notation. 
They are given by. 

Bin = {Bi^, Bin), (20) 



Cnn = I I . (21) 



where. 



Bi. = (1 - l^ewhNs + o{e) 

Bin = e(l - ^^+0 + 0{e) 

Cuu = ujmsi^-^ewhNs + oie) 
C.N = = uZj^s^ii - + oie) 

Cnn = e^ + Oi^). (22) 

After the computation of the 10 contributing one-loop diagrams, drawn in fig. ^ we find 
the analytical results presented in Appendix A which have been written in terms of the 
standard one-loop integrals, Cq, Bq,Ci2, ... etc, whose definitions can be found for instance 
in These provide the total contributions to the relevant form factors Fl and Fr that 
are related to the decay amplitude F by: 

tF = -igui^{-p,){FLPL + FRPn)viJps) (23) 
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where, 

10 10 

Fl = Y.^l\ Fn = J2F^^ (24) 

1=1 i=l 

and Pi = P3 — P2 is the ingoing Higgs boson momentum. 

The LFVHD widths can be obtained finally from these form factors by: 





mL)(|F^r + \Fr\^) - 4m,,m,„i?e(F^F*)) (25) 



Notice that since we will consider complex R matrices, the corresponding decay widths 
for the CP conjugate states, in general, can be different. We do not study here these CP 
conjugate decays and concentrate on the H rft, re, fie decays. 

Regarding the analytical results, it is worth mentioning that we have checked the finiteness 
of the total form factors Fl and Fr. When summing up all involved indices in each diagram, 
we find out, in agreement with ref. that the only divergent diagrams left are 1, 8 and 
10 and these divergences cancel among each other, providing the expected finite result. 

We next present the numerical results. The one- loop functions Cq,Bq, .. have been eval- 
uated with the Mathematica package of j^]. For the numerical estimates of the branching 
ratios we evaluate the coupling matrices of eq. EH in terms of ^, for the particular input 
values of mo and tum that are compatible with the neutrino data. That is, we get mo from 
eq. for the two scenarios A and B and for various choices of the input masses m^^^^ 
and the matrix R, as explained in sec. IHl and take iriMi ~ ^Ni- The total width has been 
evaluated with the HDECAY programme (2l|. 

The main conclusion from the results in table I is that the LFVHD rates in the SM- 
seesaw are extremely small and depend strongly on the Majorana mass scale. In all scenarios 
considered here and for large m^v, say larger than 10"^ GeV, the decay rates decrease with the 
Majorana mass, as ~ (l/m|^), and therefore they become very small for very large masses 
of the heavy neutrinos. This behaviour is easily explained by the fact that the dominant 
contributions to the form factors are proportional to ^^"^ ~ (9(m|,/m|^) and, mo being 
obtained by eq. El goes as m/j ~ ""^m^? so that the form factors scale as F ~ l/ruM and, 
consequently, the branching ratios as BR ~ l/m\.^. That is, we get decoupling of the heavy 
neutrinos. We have checked numerically this behaviour, as can be seen in fig|21 Incidentaly, 
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Ik ^k 




(1) (2) (3) 



h h h 




(7) (8) (9) (10) 

FIG. 1: One-loop diagrams for the LFVHD in the SM-seesaw model 

it should be mentioned that since we are keeping the masses of the hght neutrinos non- 
vanishing, the asymptotic value is not exactly zero but an extremely small value. 

The R matrix that appears in the relation between uid and niM can be also relevant in 
some scenarios. In particular, if the 9i angles are complex, the R-matrix elements can have 
large modulus, the derived mn can increase relevantly and, therefore, the decay rates can be 
much larger, although still very small. This can be seen clearly by comparing the first and 
second columns of table I. Notice also that, due to the leptonic mass hierarchy, the larger 
rates are always obtained for the decays involving the r-lepton. 

Therefore one can conclude that the LFVHD rates in the SM-seesaw are neghgibly small 
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Rq 


Ri 


(A) 

mjvi 2,3 = 10'^ GeV 


BRi2 
BRn 
BR23 


3 X 10'^^ 
1 X 10-39 
3 X 10^3^ 


2 X 10-32 

2 X 10-30 

3 X 10-30 


(A) 

"^iVi.2,3 = 10*^ GeV 


BR12 
BRn 
BR23 


1 X 10-"^^ 
3 X 10-^4 
1 X 10-40 


1 X 10-36 
1 X 10-34 
1 X 10-34 


(B) 

mAr^(2j = 1(5) X 10^ GeV 
ruN^ = 1 X 10^ GeV 


BR12 
BR13 

BR23 


1 X 10-3® 
1 X 10-36 
1 X 10-3' 


1 X 10-35 
3 X 10-33 
3 X 10-32 


(B) 

mAr^(2j = 1(5) X 10^ GeV 
= 1 X 10^ GeV 


BRu 
BRn 
BR23 


4 X 10-39 
1 X 10-3"^ 
1 X 10-37 


1 X 10-36 

1 X 10-34 

2 X 10-33 



TABLE I: Branching ratios of the LFVHD in the SM-seesaw, for scenarios A and B and for various 
choices of the heavy neutrino masses. The R matrix is chosen as in case and case 1 with 9i = 
92 = 03 = vr/3e''^/3. Here, BR12 = BR{H fie), BRn = BR{H re), BR23 = BR{H rp) 
and niH = 115 GeV. 



and the largest values are for the lowest Major ana mass choices and largest complex angles 
in the R matrix. These conclusions are valid for both scenarios A and B. 

In order to compare these results with the ones obtained in it has to be taken 
into account that in their computation the seesaw parameter ^ was fixed to a particular 
numerical value, so that, for large itlm-, rno scales as m^) ~ mM, in contrast to our case 
where, as we have said, rriD "^ii ■ These different assumptions on the relation between m/j 
and rriM give rise to very different decay rates. The branching ratios in their computation 
grow instead with mjvf, thus by providing very large values to the Dirac masses (i.e., the 
Yukawa couplings get strong as mjy — J^j-a/ increases) they obtained much larger values for 
the branching ratios. For a fixed neutrino mass m-Tv, the only way to reach these large values 
for the Dirac masses in our scenario is by considering extremely large complex R-matrix 
elements. As an example, in scenario A and case 1 with m^^ = 1 TeV and Q\^2,3 ~ 5z, 
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FIG. 2: Dependence of the LFVHD ratios in the SM-seesaw with mi\j (GeV) in scenario A with 
degenerate heavy neutrinos. Solid, short-dashed and long-dashed lines are for BR23, BRi^ and 
BR12 respectively. 

we obtain |m/)| ~ 200 GeV and BR23 ~ 10^^, which are closer to the predicted values in 



ref. [1^. But this is not a natural assumption in seesaw models since it requires a strong 
fine-tunning in the choice of the m£, and tum matrices. 

We also include next, for comparison, the numerical results for the LFVHD branching 
ratios in the case of Dirac neutrinos and for ttlh = 115 GeV. For the two scenarios A and B 
they are, correspondingly: 

Br{H ^ fie) = 2x 10"^^ (A), 2 x 10"^^ (B) 
Br{H ^ re) = 1 X 10"^° (A) , 1 x 10"^° {B) 
Br{H ^ r/i) = 3 X 10^^^ (A) , 2 x 10^^^ (B) 

As we can see, these ratios are negligible and much smaller than in the case of Majorana 
neutrinos. We also learn from our results than the ratios undergo a strong cancellation 
when the internal neutrinos in the loop diagrams are summed over the three generations. 
This suppression is similar to the GIM suppression mechanism of the quark sector. Notice 
that our rates disagree in many orders of magnitude with the results of ref. jl^, where a 
prediction of BR{H —>■ r/2) ~ 6 x 10~^ was presented. We believe that this disagreement 
could be due to the above mentioned strong cancellations which are not taking place in their 
case. 
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IV. LEPTON FLAVOR VIOLATING HIGGS DECAYS IN THE MSSM-SEESAW 



In this section we study the lepton flavor violating Higgs decays within the context of 
the MSSM-seesaw model. We first analyze the various sources of lepton fiavor changing 
processes in this model, and next we compute the LFVHD partial widths and branching 
ratios for the three neutral MSSM Higgs boson decays, T^H^ — > Iklm) with = ho, Hq, Aq 
and Iklm = T/i, Te, fie. We also analyze in parallel the lepton fiavor changing Ij li'-f decays 
and explore the maximum predicted rates for LFVHD, mainly for H^, rft decays, by 

requiring the BR{lj — > li'y) rates to be within the experimental allowed range. We use the 
present experimental upper bounds given, respectively, by \BR{fi — >• 67)! < 1.2 x 10^^^ ^, 
\BR{t /i7)| < 3.1 X 10-^ Q and \BR{t e-f)\ < 2.7 x 10"^ 11|. 



A. Sources of LFV interactions in the MSSM-seesaw 



In the MSSM-seesaw model there are two sources of lepton fiavor changing processes. The 
first one is induced from the non-vanishing mixing in the light neutrino sector, that is from 
the off-diagonal elements of the Umns matrix. We have already seen in the previous section 
that these elements induce, via charged currents, intergenerational interactions of the type, 
— k — Vj (and the corresponding ~ h ~ ^j) with i ^ j. These will be generated in 
the MSSM-seesaw as well, but in addition, there will appear new tree-level LFV interactions 
involving neutrinos, as for instance those with charginos and charged sleptons /, x^—li—i^j 
with i ^ j, and those with neutralinos and sneutrinos z/, — Vi — Vj with i ^ j. From the 
operational point of view, these effects are introduced by the explicit factors of the Umns 
matrix elements in the corresponding couplings when referred to the physical basis. Due to 
the extended Higgs sector of the MSSM-seesaw, there will be additional intergenerational 
couplings of the type H^Uvj with similar {UMNs)ij factors. 

All the previous intergenerational couplings can induce important contributions to LFV 
processes with external SM leptons. In the following we will refer to the radiatively induced 
fiavor changing effects from the {UMNs)ij factors in the couplings as Umns effects. For 
the case under study here of LFVHD, the couplings inducing these Umns effects are just 
W^liUj, C^liUj and H'^liUj. The size of the generated effects from the W^liUj and G^liUj 
couplings are expected to be as small as in the SM-seesaw case. The only difference comes 
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from the different couplings of the internal particles to the external ho, Hq, and Aq Higgs 
bosons as compared to those of the SM Higgs boson H, but this difference will not change 
relevantly the small size of the generated ratios. On the other hand, the generated effects 
from the H^liUj couplings are new with respect to the SM-seesaw case. We have numerically 
estimated the size of the contributions to the LFVHD branching ratios from all the one-loop 
diagrams with internal H"^, and G^, and we have found that they are indeed always 
negligibly small. We will in consequence ignore all these contributions in the following. 

The second source of LFV processes in the MSSM-seesaw is genuine in SUSY models 
and comes from the misalignment between the rotations leading to the mass eigenstate 
basis of sleptons relatively to that of leptons. This misalignment is radiatively generated 
in the SUSY-seesaw models from the Yukawa couplings of the Majorana neutrinos and can 
be sizeable in both the charged slepton and sneutrino sectors . Once one rotates the so- 
generated flavor non-diagonal charged slepton and sneutrino mass matrices to the physical 
diagonal ones, some intergenerational couplings involving SUSY particles are generated. For 
the case of LFVHD, the involved couplings are x^Uvj and j<^lilj with i ^ j. These generated 
effects from lepton-slepton misalignment will be called in the following misalignment effects. 
All these couplings will induce via loops of SUSY particles relevant contributions to the 
LFVHD rates as will be shown later on. 

The LFV effects from misalignment are usually implemented in the seesaw models in the 
language of the Renormalization Group Equations (RGE) and can be summarized as follows. 
One starts with universal soft-SUSY-breaking parameters at the large energies Mx » tum 
given by, 

{mi)l = Ml5i^, {m^)l = M^S^j, (m^)J = M^S^j, (A^),,- = Ao{Yi)ij, (A,),,- = Ao{Y,U26) 

where, Mq and Aq are the universal soft-slepton mass and soft-trilinear coupling respectively. 
Yi is the Yukawa coupling matrix of the charged leptons, which is flavor diagonal in the basis 
chosen here, {Yi)ij = Yi^6ij with Yi. = ^ and vi = v cos (3. Y^ is the Yukawa coupling matrix 
of the neutrinos and is given by, iYy)ij = , with f 2 = f sin f3. 

The effect of the RGE-running from Mx down to tjim on the off-diagonal soft parameters 
of the charged slepton sector, to one-loop and in the leading-log approximation, is then 
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described by, 

(Am|),, = --L(3M2 + A2)(F;LrJ),, 

(Ami),, = 0;Lki = log ( ^) 5^. (27) 

Notice that our notation is slightly different from the usual one in that i — > Yj. For 
all the numerical estimates in this work we use Mx = 2 x 10^^ GeV and, for simplicity, 
we will assume Aq = 0. Thus, we only consider flavor changing in the LL sector, which 
is known to be a very good approximation within the context of the seesaw model. In 
fact, we have checked that even for values as large as Aq = Mq = ITeV the size of the 
corresponding flavor changing dimensionless parameters in the LR sector are always smaller 
than the LL ones in more than three orders of magnitude. Notice also that, in addition to the 
previously mentioned Umns factors, which, as we have said, are not relevant for the present 
computation of the LFVHD rates, there is an extra dependence on Umns via the Y^, couplings 
or, equivalently, via the rriD matrix as given in eq. EH Thus, even if we fixed Umns to the 
unit matrix, there could be still flavor changing effects from misalignment via the R matrix. 
Conversely, if we fixed R to the unit matrix, there could be still flavor changing effects 
from misalignment via the Umns matrix in m^,. Finally, notice that the effect of neutrino 
Yukawa couplings, via RGE-running from Mx down to tum, on the diagonal entries of the 
squared-mass matrices are small, and have been neglected in this work. On the other hand, 
as we have previously said, we have included the effects of the RGE-running from mw up 
to rriN on the neutrino masses m^- {i = 1, 2, 3) and the Umns matrix elements, and we have 
used the corresponding corrected parameters in eq. ^Ito evaluate the Yukawa couplings Y^. 

In order to estimate the size of the misalignment effects, we first study in detail the 
dependence of the flavor changing dimensionless parameters, in the LL sector, defined as 
\ij = ^—j^-^, as a function of the seesaw parameters. We show in figs. Q and (H)) 
the predictions of the Aj^ parameter as a function of the seesaw parameters, and in some 
selected examples within the scenarios described above. We show in fig. (jHt^) the dependence 
of \Xij\ with rriN, for scenario A and real R. Notice that in this case, the Xij turn out to be 
independent on R, so this figure applies both to case and to case 1 with 6i real. As can 
be seen, the three \Xij\ grow with m^v and the largest one, which is IA23I, reaches values up 
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FIG. 3: Dependence of Xij with the seesaw parameters. Sohd, long-dashed and short-dashed hnes 
are for A23, A13 and A12 respectively. (3a) Upper left panel: Dependence with m-Ar in scenario A for 
real R. (3b) Upper right panel: Dependence with the a parameter in scenario A with mj\[ = 10^^ 
GeV and complex R, for case 2 with a = b = c. (3c) Lower left panel: Dependence with real 9i 
in scenario B, case 1, for {1711^^,1111^1^,711^.^) = (10^,2 x 10^,10^^) GeV and 62 = 9z = 0.(3d) Lower 
right panel: Same as in (3c) but for Umns = 1- In all plots here, tan/3 = 35. 

to 2.4 X 10"'^ for m^r = 10^^ GeV. The size of IA12I can reach values up to 3.8 x 10^^ for 
rriN = 10^^ GeV, and correspondingly for |Ai3| ~ |Ai2|. The size of the Xij in this scenario 
with degenerate heavy neutrinos can obviously be increased if R is assumed instead to be 
complex. In this case, the Xij are in general complex numbers. As an example, we show in 
fig- (Eb) the dependence of \Xij\, in scenario A with rriN = 10^^ GeV, on the parameter a of 
case 2, for a = 6 = c . The size of \Xij \ increases clearly with a and, for the studied range, IA12I 
can be as large as 2 x 10^^. In particular, for values of \abc\ ^ 10^^, that generate succesfuU 
baryogenesis [q], the generated \ Xij\ are still large, namely, IA12I = 4 x 10^^. Notice also that 
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FIG. 4: Dependence of \Xij\ with the seesaw parameters for complex 6i. Solid, long-dashed and 
short-dashed lines (the two later being undistinguishible in the plots) are for IA23I, |Ai3| and |Ai2| 
respectively. Left panels: Dependence with \9i\, \92\ and |03| respectively in scenario B, case 
1. We take Arg{9i) = 7r/4, correspondingly, and the rest of angles are set to zero. Right panels: 

Dependence with Arg{6i), Arg{92) and Arg{9^) respectively in scenario B, case 1. We take \9i\ = vr, 
correspondingly, and the rest of angles are set to zero. In all plots here, {niN-i^, 171^2, n^Na) = 
(10^ 2 X 10^ 10^^) GeV and tan/3 = 35. 
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the relative size of the different Xij changes with respect to fig. ^.). The main problem 
with this scenario A and R complex, with the parametrization of case 2, is that for most 
choices of the relevant parameters, m^, a, b and c, producing succesfuU baryogenesis, the 
size of the generated A12 is too large and leads to large yU — > 67 

the present experimental bound of BR{fi 67) < 1.2 x 10~^^ [9,]. We have also studied 
the scenario A with complex R given by the parametrization of case 1, and have found that 
for most choices of the 9i complex angles the three A12, A13 and A23 have comparable sizes 
and therefore, the required large A23 values leading to large H ^ rfl ratios imply again 
experimentally unallowed /i — >^ 67 ratios. Therefore, we will not consider this scenario A 
with complex R anymore in the following. 

The case of hierarchical neutrinos is shown in figs. (jSt) and (jSll) for real R and in fig. (@)) 
for complex R with the parametrization of case 1. In fig. (jSfc) it is plotted the dependence 
with real 61 for scenario B with {rriN-^^, m^rj, m^rg) = (10^, 2 x 10^, 10^^) GeV and 02 = 03 = 0. 
We see that A23 is the largest one and reaches negative values up to —1.96 x 10^^ at the 
points 01 = 0, ±7r, precisely where A12 and A13 vanish. These points are therefore particularly 
interesting for the LFVHD rates since they will lead to the largest ratios for H rp, while 
keeping fi —>■ 67, r ^ 67 extremely small and, as we will see later, t ^ fi'-f still compatible 
with data. Notice that the point 6'i = in fig. (jSh) corresponds to the reference case 
with R = 1 and, therefore, represents the situation where the Umns niatrix is the only 
origin for flavor changing. This means that our experimental input for the Umns niatrix 
generates by itself sizeable rates for lepton flavor violating decays involving the second and 
third generations. The alternative situation where just the R matrix is generating the lepton 
flavor violating decays is illustrated in fig. (jSli). Here we show the Xij dependence with 0i 
for 6'2 = 6'3 = and Umns = 1- We see that IA23I reaches values up to about 7.8 x 10~^ 
whereas A12 and A13 vanish for all 0i. By comparing these two situations we can infere that, 
for the case of real R with 6*1 7^ and 02 = 03 = 0, the induced misalignment effect between 
the second and third generations from the experimental Umns is relevant and can be even 
larger than the effect from R. 

We have also studied the alternative choices for real R with ^^2 7^ 0, = ^3 = and with 
^^3 0, 01 = 02 = 0, although the corresponding plots are not shown here for brevity. We 
find a A23 dependence on 02 very similar to that on 6*1, with maximum negative A23 values 
at 02 = 0, ±7r of —1.96 x 10^^. In contrast, A12 and A13 now take very small values whose 



21 



maximum are 2.4 x 10~^. Regarding the dependence with ^3 a different situation is found, 
where the three A12, A13 and A23 are approximately constant with ^3 and take the values, 
A23 = —1-96 X 10"^, and A13 ~ A12 = 2.4 x 10"^, respectively. 

The case of hierarchical neutrinos with complex R produces, in most cases, complex Ajj 
values and their moduli are in general larger than in the case of real R, as can be clearly 
seen in fig. (jl)). Our particular choice for the heavy neutrino masses of {rriN^, 171^2, ^n^) = 
(10^, 2 X 10^, 10^^) GeV, where the two lightest neutrinos have similar masses and well below 
the mass of the heaviest one, produces the specific pattern shown in these plots, where the 
dependence of IA23I on 6*1, for 6^2 = 6*3 = 0, and on 6*2, for = ^^3 = are very similar, and 
IA23I can reach very large values for a large region of the {\6i\,Arg{6i)), i = 1,2 parameter 
space. For instance, for fixed Arg{9i) = 7r/4, and |6'i| up to n we find IA23I values up to 
0.46 and similarly for 62- Larger values of Arg{6i), i = 1,2, produce even larger IA23I and it 
reaches its maximum at Arg{6i) = tt/2. In contrast, A12 and A13 reach much smaller values 
with complex 62 than with complex 61, being IA12I — IA13I < 5 x 10^^ for \62\ < vr. On the 
other hand, they depend strongly with complex 61 and IA12I can reach too large values, up 
to O{10^^), in clear conflict with the allowed values by /i — ^ 67 data. Finally, the behaviour 
with complex 6*3 is very similar to the real case, with the three IA12I, IA13I and IA23I being 
nearly constant with 6^. Their values are IA23I = 1.96 x 10~^ and IA12I — IA13I = 2.4 x 10~^, 
respectively. 

Regarding the values of the hierarchical heavy neutrino masses, we have also tried other 
choices in the range lO^GeV < rriN. < 10^'^GeV and found that, under the favored assump- 
tion by baryogenesis of close rriNi and 171^2 and much lighter than 171^3, it is the value of 
this later what matters for BR{H rfl), leading to larger A23 values for larger m^.^ and, 
therefore, our choice of {771^-^,771^^,171^,^^) = (10^,2 x 10^,10^'') GeV seems to be appropri- 
ate. We have checked that alternative choices where m^Va and m^vg are close and much 
heavier than ttln-^ lead to a similar situation for 61 ^ 0, 62 = 63 = than the previous 
case, but it gives larger preditions for A12 in conflict with fi e'-f data. For the follow- 
ing predictions of decay rates in the case of hierarchical neutrinos we will fix the values to 
(mjVi,m^2,^iV3) = (10^2 X 10^ lO^^) GeV. 

In summary, the case of hierarchical heavy neutrinos leads to larger Xij values than the 
degenerate case and, in consequence, larger LFVHD rates. Futhermore, in order to get 
the largest possible BR{H — » rfx) rates while keeping all BR{lj li'-j) rates within the 



22 



experimental allowed regions, the most favorable case of all the studied ones, and for our 
choice of neutrino masses, is the one with complex ^2 7^ 0, and = = 0. 

The next stept is the diagonalization of the charged slepton and sneutrino mass matrices 
leading to the mass eigenvalues and mass eigenstates at the electroweak energy scale. We 
start with the non-diagonal charged slepton and sneutrino squared-mass matrices that are 
obtained after the running from Mx to Mw and once the charged leptons and neutrinos 
have been rotated to the physical basis. For the charged sector, this matrix is referred to 
the (e^, cr, p^L, flfi, tl, tr) basis and can be written as follows. 
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The soft SUSY breaking masses and trilinear couplings above, ^, mi ^, ^ and Ai, refer 
to their corresponding values at the electroweak scale. We got them by solving numerically 
the RGE with the program mSUSPECT [2^ and by imposing the universality conditions at 
Mx = 2 X 10^^ GeV for the sfermion sector, eq. EHl together with the corresponding ones 
for the gaugino and Higgs boson sectors, Mi(Mx) = M2{Mx) = M1/2 and Mh^{Mx) = 
Mh^{Mx) = Mq respectively. For the gaugino sector, this implies the well known relation at 
low energies, Mi(Mvy) = ^{ta.n 9wf M2{Mw) ■ The value of the supersymmetric fi parameter 
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is extracted as usual from the electroweak breaking condition. We choose in all this paper, 
/i > 0, and do not expect relevant differences for fi < 0. 

After diagonalization of the M? matrix one gets the physical slepton masses and the 

six mass eigenstates (/i, ,^6)= ^ which are related to the previous weak eigenstates 

(eL,....fji)= I' by the corresponding 6x6 rotation matrix, /' = R^'''l. 

Regarding the sneutrino sector one proceeds similarly to the charged slepton sector, 
but now the diagonalization process is simpler because of the involved seesaw matrix ^ 
which gives rise naturally to a suppresion of the RH sneutrino components in the relevant 
mass eigenstates. This can be easily illustrated in the one generation case, but for three 
generations one arises to similar conclusions. The sneutrino mass terms of the MSSM-seesaw 
.„de.canbew«e„in*e„nege„e.a«o„ca.eHa3. 
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with. 
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? + m^) — im| cos2/5 m a {A^ — fi cot jS ± niM] 



(31) 



mo (Ay - /i cot /? ± itim) rri'^^j + m|, + ± 2BMmM 

Notice that now there are several mass scales involved, the soft SUSY-breaking param- 
eters, m^, m^, Bm and A;^, the Dirac mass m^), the /x-mass parameter, the Z boson 
mass mz and the Majorana neutrino mass itim- Our basic assumption in all this paper 
is that mM is much heavier than the other mass scales involved (obviously, except Mx), 
niM » ^D,fnz,fi',^L^^M^^i^^BM- The size of Bm has been discussed in the litera- 
ture 23] and seems more controversial. For simplicity, we shall assume here this is also 
smaller than rriM- In this situation, the diagonalization of the previous sneutrino squared 
mass matrix is simpler and leads to four mass eigenstates, two of which are light, C,{, C,2 and 
two heavy, C,i, C,2- In the leading orders of the series expansion in powers of ^ the mass 
eigenstates and their corresponding mass eigenvalues are given by. 
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(32) 



Here we can see that the heavy states 2 ^^^^ couple very weakly to the rest of particles 
of the MSSM via their i>L component, which is highly suppresed by the small factor ^ and, 
therefore, it is a good approximation to ignore them and keep just the light states ^1 2) which 
are made mainly of i>L and its complex conjugate z>£. Now, by working in this simplified 
basis but applied to the three generations case, which we write for short i>'^ {a — 1, 2, 3), the 
relevant 3x3 sneutrino squared mass matrix can be written as follows. 
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where m?^ and (Am?)ij are the same as in the previous charged slepton squared mass 
matrix. After diagonalization of the M? matrix one gets the relevant physical sneutrino 
masses and eigenstates, (/3 = 1, 2, 3) which are related to the previous states z/^ by the 
corresponding 3x3 rotation matrix, u' = R^'^^u. 

To end up this subsection, we summarize all the interaction terms that are relevant for 
the computation of the LFVHD rates. We present these interactions in the physical mass 
eigenstate basis and will perform all the computations in this basis. It implies diagonaliza- 
tion in all the involved SUSY sectors, charged sleptons, sneutrinos, charginos, neutralinos 
and Higgs bosons. The SUSY-electroweak interaction terms among charginos, leptons and 
sneutrinos and among neutralinos, leptons and sleptons, that are the responsible for the 
LFVHD are as follows. 
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where the coupling factors A^2xj, A^^j, B^H^^, and B^^^ are given in Appendix B. 
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(34) 
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The other interaction terms that enter in the computation of the LFVHD rates are the 
Higgs-lepton-lepton, Higgs-sneutrino-sneutrino, Higgs-slepton-slepton, Higgs-chargino- 
chargino, and Higgs-neutrahno-neutrahno interactions, reading: 



HxSaS/s 



-gHJ 

-iH. 



^xXaXt 



where the couphng factors ^^g. S^ni^ Qh^u^v^. 9hJJp^ ^uj, Wmj^ D'^ab and D)^^^ are col- 
lected in Appendix B. 
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B. LFVHD rates in the MSSM- seesaw 



As we have said, the contributions to the LFVHD rates in the MSSM-seesaw come from 
various sectors. The contributions from the charged Higgs sector and from the SM sector 
(i.e, if^, and G^) are very small and will not be included here. The main contributions 
come from the genuine SUSY sector, concretely, from the one-loop diagrams with charginos, 
neutralinos, sleptons and sneutrinos shown in fig. |S| 

The contributions of these one-loop diagrams to the form factors are given by. 



i=l 



(36) 



i=l 



where the analytical results for and , F'^'^, i = 1..8, are collected in Appendix B. 

The partial widths for the H^, A° Iklm decays are then finally computed by inserting 
these form factors correspondingly into eq. ESI We show in figs. Elto^lthe numerical results 
of the branching ratios for the LFVHD in the MSSM. The total MSSM Higgs boson widths 
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211 . We have shown in the plots just 



have been computed with the HDECAY programme 
the dominant channels, which are — > t/2, and some comments will be added on the other 
channels. Similarly, for the comparison with the leptonic radiative decays, Ij k'y, we will 
show in the plots the most relevant one, which is /i — > 67 or r ^ fij, depending on the case. 
For the numerical estimates of the Ij li'y branching ratios we use the exact analytical 
formulas of ref. These are expressed in the mass eigenstate basis as well, and contain 
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(5) (6) (7) (8) 

FIG. 5: One-loop diagrams for the LFVHD in the MSSM-seesaw model 

all contributing one-loop diagrams. For the involved x~^^ couplings we use again 

the expressions of Appendix B. 

The results of the branching ratios for the LFVHD, in the r/i channel, as a function of 
the Majorana mass, m^v, in scenario A with degenerate heavy neutrinos and real R, are 
illustrated in fig. El for several tan/5 values, tan (3 = 3,10,30,50. The explored range in 
is from 10^ GeV up to 10^^ GeV which is favorable for baryogenesis. We also show in 
this figure, the corresponding predicted rates for the most relevant lepton decay, which in 
this case is /i — 67, and include its upper experimental bound. We have checked that the 
other channels are well within their experimental allowed range. From this figure we first 
see that the branching ratios for the light Higgs boson are smaller than the heavy Higgs 
ones in about two orders of magnitude. The ratios of Hq and Aq are very similar in all 
the plots and, for this scenario, they can reach values up to just 2.2 x 10~^° in the region 
of high tan/3 and high m^r. Besides, the rates for /i — > 67 decays are below the upper 
experimental bound for all explored tan/5 and m^r values. From these plots we also see 
clearly the high sensitivity to tan (3 of the LFVHD rates for all Higgs bosons which, at large 
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FIG. 6: Dependence of BR{Hx tJI) with ttin (GeV) in scenario A with degenerate heavy 
neutrinos and real R, for several values of tan/3. (6a) Upper- left panel, Hx = Hq, (6c) lower-left 
panel, Hx = Aq and (6d) lower-right panel, Hx = h^. (6b) Upper-right panel: Dependence of 
BR{iJ, —>■ 67) with mjv for several values of tan/3. In all plots, the solid, dashed, dashed-dotted 
and dotted lines are the preditions for tan /? = 3, 10, 30 and 50, respectively. The horizontal line in 
(6b) is the upper experimental bound on BR{iJ, ej) . The other input parameters are, Mq = 400 
GeV and M1/2 = 300 GeV. 

tan/3, scale roughly as (tan/3)'^, in comparison with the lepton decay rates which scale as 
(tan/3)^. The dependence of both rates on mjv is that induced from the Xij dependence, 
and corresponds approximately to what is expected from the mass insertion approximation, 
where BR{Hx Ijh), BR{lj lij) oc \Xij\'^ oc |mivlog(mjv)p. 

In what regards to the relative importance of the various SUSY sectors to the LFVHD 
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rates, we have found that these are dominated by the chargino contributions, that is from the 
loop diagrams (1),(2),(3) and (4) in fig. El For instance, for niN = 10^'' GeV , Mq = 400 GeV, 
and Ml/2 = 300 GeV, we have found the following ratios between the chargino and neutralino 
contributions to the form factors: {Pf /Ff\ = 6.1, 6, 7.3, 23.1 for tan/3 = 3, 10, 30, 50 
respectively, where we have used a simplified notation, = Pj^l^^^ + Pj^lj^ + Pj^lj^ + F^]^^, 
= f]^^^^ + Ff^^^ + i^i'l'l/Q + Fl'^Hq- Similar ratios are found for the corresponding 

Fn form factors. The relative ratio found of F^/ F^ ~ 17 is nicely explained by the m^/m^ 
ratio. For the lightest Higgs boson, we find \F^ / F^ \ = 1.5, 1.4, 1.7, 4 correspondingly. 

Concerning to the comparative size of the contributions from the various chargino loop 
diagrams we have found that, at large tan/3, the external leg corrections are clearly the 
dominant ones. Concretely, for \iFi%^+ Fi%^)/Ff\, \Fi%jFf\ and \Fi%jFf\, we get 
the respective percentages, 60.6%, 39.3% and 0.1%, for tan/3 = 10 and 93.8%, 6.2%, 0% for 
tan (3 = 50. 

The branching ratios for the Higgs boson decays into re and fie are much smaller than 
the r/2 ones, as expected, and we do not show plots for them. For instance, for rriN = 10^^ 
GeV, and tan/3 = 50 we find BR{H'^'''> Tjj,)/BR{H^''^ re) = 3.9 x 10^ and BR{H^'''^ 
Tfi)/BR{H'-''^ fie) = 1.3 X 10^ for the three Higgs bosons. 

All the previous results are for fixed Mq = 400 GeV and M1/2 = 300 GeV. The dependence 
with Mq and M1/2 will be discussed later on within the context of hierarchical neutrinos. 

In summary, the LFVHD rates for degenerate heavy neutrinos are very small, at most 
2.2 X 10~^°, for the explored range of the seesaw parameters and tan/3. Obviously, larger 
values of these LFVHD ratios could be got for larger tan /? values, but we have not considered 
them here. 

We next present the results for hierarchical neutrinos, scenario B, and use the 
parametrization of case 1. The results for real and complex R and for the mass hierar- 
chy {mN^.mN^.rriN.^) = (10^,2 x 10^,10-'^^) GeV, are shown in figs, d and |H1 respectively. 
From these figures we first confirm that the LFVHD and lepton decay rates are larger in 
this case than in the degenerate heavy neutrinos one. However, we will get restrictions 
on the maximum allowed Higgs decay rates coming from the experimental lepton decay 
bounds. For instance, the case of real 9i, that is illustrated in figures ((7^) and ((7|d) for 
tan/3 = 50, Mq = 400 GeV and M1/2 = 300 GeV, shows that compatibility with fi ^ e'j 
data occurs only in the very narrow deeps at around 61 = 0, 1.9 and vr. The presence of 
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FIG. 7: (7a) Left panel: Dependence of BR{Hx — > rji) with ^i. Solid, dashed and dashed-dotted 
lines (the two later undistinguishible here) correspond to = {ho, Hq, Aq) respectively. (7b) 
Right panel: Dependence of BR{fi 67) with 61. The horizontal dashed-dotted line is the upper 
experimental bound. Both panels are in scenario B, case 1 for real 9i 7^ 0, (m-TVi , ?T^Ar2 ) ""^Afa ) = 
(10^,2 X 10^10l4) GeV, 62 = 63 = 0, tan/3 = 50, Mq = 400 GeV, and M1/2 = 300 GeV. (7c), 
lower left panel, and (7d), lower right panel, are as in (7a) and (7b) respectively, but for ^2 7^ 
and 6*1 = 6*3 = 

these narrow regions where the /i — 67 rates are drastically suppresed were already pointed 
out in ref. and correspond clearly to the minima of IA12I in fig. (j2t). Notice that it is 
precisely at the points 6'i = 0, vr where the BR{Hq, Aq rft) rates reach their maximum 
values, although these are not large, just about 1.3 x 10^*^. Notice also, that these maxima 
correspond clearly to the maxima of IA23I in fig. (jHt). We have checked that for lower tan/? 
values, the allowed regions in 9i widen and are placed at the same points, but the corre- 
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FIG. 8: (8a) Upper-left panel: Dependence of BR{H^ r/x) with Arg{e2) for I6I2I = tt. (8b) 
Upper-right panel: Dependence of BR{t ^7) with Arg{62) for |02| = t^- (8c) Lower-left panel: 
Dependence of BR[Hx Tp) with |02| for Arg[92) = '7r/4. (8d) Lower-right panel: Dependence 
of BR[t fij) with|02| for Arg{92) = 'tt/4:. All figures are in scenario B, case 1, for complex 82 
and 9i = 9-^ = 0. The rest of parameters are fixed to: {mj^i^^mj^i^^mN^i) = (10^, 2 x 10^, 10^^) GeV, 
tan 13 = 50, Mq = 400 GeV and M1/2 = 300 GeV. SoHd, dashed and dashed-dotted (the two later 
undistinguishible here) lines in the left panels correspond to = {ho, Hq, Aq) respectively. The 
horizontal line in the right panels is the experimental upper bound on r ^ fi^. 

sponding maximum values of the LFVHD rates get considerably reduced. The alternative 
case of real 62 7^ 0, with 6'i = 6*3 = is illustrated in figs. ^) and (I7|i). We see that 
the behaviour of BR{Hx rfl) with 62 is very similar to that with 61 of fig. (|7^) and the 
maximum values of about 1.3 x 10~^ are now placed at 62 = 0, vr. BR{fi 67) also reaches 
its maximum at 62 = 0, tt, but it is still well below the experimental bound. In particular, 
for tan/3 = 50, Mq = 400 GeV and M1/2 = 300 GeV this maximum value is 3 x 10"^^ 
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Notice, that the behaviour with 62 is explained once again in terms of the corresponding Xij 
behaviour. Regarding the dependence with 6^, not shown in the plots, a different situation 
is found, where BR{Hx — » r/2) is approximately constant, and for the heavy Higgs bosons it 
is around 1.3 x 10^^. BR{fi 67), BR{t fi'-f) and BR{t —>■ 67) are also approximately 
constant with 6^. In addition, we have checked that these three leptonic constant decay 
rates are within the experimental allowed range. In summary, for real R we find that the 
maximum allowed LFVHD rates are at or below 1.3 x 10~^. 

The case of complex R is certainly more promissing. The examples shown in figs. (jHK) 
and (jHfc) are for the most favourable case, among the ones studied here, of complex 7^ 
with 61 = 63 = 0. It shows that considerably larger BR{Hx — > tJi) rates than in the real 
R case are found. For the explored 62 values in these plots, the Higgs rates grow with both 
1^21 and Arg{62) and, for the selected values of the parameters in this figure, they reach 
values up to around 5 x 10~^. We have checked that the predicted rates for BR{t 67) 
are well below the experimental upper bound and that the fi ^ e'y decay is, in this case, 
less restrictive than the r ^ fi'j decay. Notice that the smallness of the fi —>■ ej and r — > 67 
decay rates, in the case under study of 7^ 0, is not maintained if our hypothesis on 
6'i3 = is changed. For instance, for ^13 = 5", which is also allowed by neutrino data, we 
get BR{n 67) ~ 2.4 x 10~^, for 62 = vre*^, well above the experimental upper bound. 
This is why we keep ^13 = in all this work. Therefore, in this case of complex ^^2 7^ 
with ^13 = 0, the relevant lepton decay is r yU7 which is illustrated in figs. (jHb) and 
dHli) together with its experimental bound. We see that the allowed region by r ^ /i7 data 
of the {\92\, Arg{92)) parameter space implies a reduction in the Higgs rates, leading to a 
maximum allowed value of just 5 x 10^^. 

These results are for fixed values of {171^^^, 771^2, ^Ng) = (10*^,2 x 10®, 10^'') GeV, tan/3 = 
50, Mo = 400 GeV and M1/2 = 300 GeV. We have found that other choices of the soft SUSY 
breaking mass parameters, Mq and M1/2, are more efficient in order to get larger maximum 
allowed Higgs ratios. For instance, for Mq = M1/2 = 1200 GeV, we find maximum allowed 
values of around 5.6 x 10^®. The reason for this improvement is the different behaviour with 
these parameters of the LFVHD and the lepton decay rates, which will be studied in more 
detail next. 

In fig. iniwe show the dependence of BR{H^ —>■ tJi) and BR{t — ^ fi'-f) with Mq and M1/2 
for hierarchical neutrinos with [171^1,1^^2, f^Na) = (10^,2 x 10®, 10^^) GeV and fixed values 
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BR{Ho Tjl) BR{t fi-f) 




FIG. 9: Dependence with Mq (GeV) and M1/2 (GeV) for scenario B with {mNijmNi^fnNs) = 

(10^2 X 10^10^'') GeV, 6'i = 6'3 = and tan/3 = 50. (9a) Upper left paneh Behaviour of 

BR{Ho Tp) with Mo(GeF) for M1/2 = 300 GeV and 62 = 7re°•^^ (9b) Upper right panel: 

Same as (9a) but for BR{t fi-f). (9c) Medium left panel: Behaviour of BR{Hq — *• rp,) with 

Mi/2{GeV) for Mq = 400 GeV and 62 = 7re°-^*. (9d) Medium right paneh Same as (9c) but for 

BR{t n'j). (9e) Lower left panel: Behaviour of BR{Ho tJI) with Mq = Mi/2(GeF) for 

62 = 7re°-^*. (9f) Lower right panel: Same as (9e) but for BR{t — /X7). The horizontal line in the 

right panels is the upper experimental bound on BR{t jjL'y). 
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of tan (3 = 50, 62 7^ 0, and 61 = 63 = 0. We see clearly in these plots the different behaviour 
of these two observables with the soft SUSY breaking mass parameters. Figs. (jHK) and (j^b) 
show a milder depencence of BR{H^ rpi) on Mq and M1/2 than that of BR{t yU7) in 
figs. (jHb) and (jntl) respectively. This implies, that for large enough values of Mq or M1/2 
or both the BR{t fi'y) rates get considerably suppresed, due to the decoupling of the 
heavy SUSY particles in the loops, and enter into the allowed region by data, whereas the 
BR[H^ Tp,) rates are not much reduced. In fact, we see in figs. (jHt) and (jHF) that for 
the choice Mq = M1/2 the tau decay ratio crosses down the upper experimental bound at 
around Mq = 1200 GeV whereas the Higgs decay ratio is still quite large ~ 6 x 10~^ in 
the high Mq region, around Mq ~ 2000 GeV. This behaviour with the soft SUSY breaking 
parameters is a clear indication that the heavy SUSY particles in the loops do not decouple 
in the LFVHD, in much the same way as it has been shown to happen in the case of Higgs 

non-decoupling of the SUSY 

particles in the LFVHD can also be reformulated as non-decoupling in the effective H^^Wfj, 
couplings and these in turn can induce large contributions to other LFV processes that are 
mediated by Higgs exchange as, for instance, r [2^. However, we have checked 

that for the explored values in this work of Mq, M1/2, tan/3, R and mi^. that lead to the 
anounced LFVHD ratios of about 6 x 10~^, the corresponding BR{t ^fi^) rates are below 
the present experimental upper bound of 2 x 10~^ |27 |. 

Finally, in order to show more clearly the non-decoupling behaviour of the SUSY particles 
in the contributing loops to the LFVHD we consider, instead of mSUGRA, a simpler and 
more generic MSSM scenario, with the Xij being free parameters, which we now fix to some 
particular values, concretely A23 = —0.4, and A12 = A13 = 0. For simplicity, we also assume 
a common SUSY mass at the electroweak scale, Msusy = ^ii = ^ei = ^0 = A* and 
choose M2 = Ml = I tan^ 6'vi/M2. This particular value of A23 corresponds roughly to 
the predicted A23 in the MSSM-seesaw with the parameters set in fig. (jHt). Finally, the 
BR{Hq —>■ Tfi) is shown in fig. El as a function of this common Msusy scale, for tan/3 = 50 
and rriHo = 340 GeV. We see clearly that for large Msusy the branching ratio approaches 
to a constant non- vanishing value, which for these input parameter values is of about 10"^, 
and therefore the charginos, neutralinos, charged sleptons and sneutrinos do no decouple in 
this observable. Another way of seing this explicitely, is by the analytical computation of 
this observable in the large SUSY masses limit. We have performed this computation in the 
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simplest case where all the SUSY masses are equal and got the following asymptotic limits 
for the dominant form factor F/,, in the regime of small A23 and large tan/3: 



(HO) 



a 



Att sin 9w ^'^rn^ 
a m-r 



A23(tan/5)2 

(1 -3tan^ ew)X23{tSiia(3y 



(37) 
(38) 



sin 6w 24miy 

From these simple expressions we can estimate quite easily the LFVHD ratios. For instance, 
for the parameters chosen in fig. we get BR{H^ — > rfx) ~ 3 x 10^® in reasonable 
agreement with our numerical result in this figure and in fig. (jSh). Notice also that this 
asymptotic result agrees with the result from the effective lagrangian approach in ref. jlfil ]. 
Finally, it is worth mentioning that this non-decoupling behaviour is in contrast with the 
behaviour of BR{ji — > 67), which scales as {My/ /MsusyY, and explains the comparatively 
large LFVHD rates found here. 
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FIG. 10: Behaviour of Hq ^ r^u in a generic MSSM scenario as a function of the common SUSY 
mass, MsusriGe^y) = i = it^e i ~ ^0 = ^- The gaugino soft masses are set to M2 = 2/3/i and 
Ml = I tan^ ewM2. Here we fix A23 = -0.4, A12 = A13 = 0, tan (3 = 50 and m^o = 340 GeV. 



V. CONCLUSIONS 

In this paper we have studied in full detail the lepton flavor violating Higgs boson decays 
that are produced if the neutrinos get their masses via the seesaw mechanism. We have 
considered the two most popular seesaw models with three generations, the SM-seesaw and 
the MSSM-seesaw. Within the SM-seesaw we have found extremely small branching ratios 
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which are explained in terms of the decouphng behaviour of the heavy Majorana neutrinos 
and the smallness of the hght neutrino masses. In the MSSM-seesaw we find, in contrast, 
branching ratios that are many orders of magnitude larger. The larger ratios found are 
for Hq Tfl and Aq r/2 decays with similar rates. After exploring the dependence of 
the Hq — s> r/i decay rates with all the involved parameters of the MSSM-seesaw, and by 
requiring compatibility with data of the correlated predictions for /i 67, r ^ 67 and 
T ^ fij decays, we find that BR{Hq —>■ rfi) as large as 10~^, for hierarchical neutrinos and 
large Msusy can be reached. These ratios are mostly sensitive to tan f3, the heaviest neutrino 
mass mATg and the complex angle 62, which have been taken in the range 3 < tan/5 < 50, 
lO^GeV < TTiN, < lO^^GeV and (l^al, ^^^(^2)) < (3.5,1) respectively. The largest allowed 
ratios found in this work of about 10^^ are for tan/3 = 50, m^s = 10^'^ GeV, large Msusy 
in the TeV range and for our choice of 62 = TTe^'^\ di = 63 = 0, but a more refined analysis 
of the full parameter space could lead to even larger rates. In particular, it is obvious that 
larger tan 3 values will enhance considerably the rates and lead to Higgs ratios closer to the 

n n 

future experimental reach of 10 ^ at LHC [2^ and e~^e and jj'^jj colliders [2^, but we 
have not tried this because it would require to perform a resummation of the large tan (3 
contributions that is beyond the scope of this work. 
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APPENDIX A 



We present here the analytical results for the form factors in the SM-seesaw, in the 
Feynman 't Hooft gauge. 

-.2 
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where Co,ii,i2 = Co,ii,i2(mf^, m|^, m^., m^, m^) and Co = Co{mf^,mjj, m^., m^, m^). 
where Co,ii,i2 = Co,ii,i2(mf^, m|^, m^., m^, m^). 
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where So,i = So,i(mf^, m^., m^). 

C'o(P2>Pi>^i>^2>^3) = Bo{pj,rnl,ml) +rnlCo{pl,pj,mj,rnl,ml) 

In all the previous formulas, summation over all indices are understood. These run 
i,j — 1, ..6 for neutrinos, and A;, m = 1, ..3, for charged leptons. 
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APPENDIX B 



1. Couplings in the MSSM-seesaw 



We 



present here the couphng factors entering in the MSSM-seesaw formulas. 
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^ — sin 



and a 



(x) 



COS a 




^ cos a ^ 



0"^' = sin a , for x = respectively, 
y u y \i cos /5y 

The matrices that rotate to the mass eigenstate basis are: U and V for charginos, for 
neutrahnos, i?*-'-' for charged sleptons and R^'^^ for sneutrinos. f/, V and are taken from 
ref. ^ai = ^ai cos 9w + Na2 sin6'vF, A^a2 = -^ai sin6'iy + A^a2 cos9w, and and i?^'') 
are computed here by the diagonahzation procedure presented in section IV. The various 
indices in the previous formulas run as follows: i,j = 1,2 for charginos, a,b = 1,..,4 for 
neutralinos, = 1, ..,6 for charged sleptons, = 1, ..,3 for sneutrinos, and I = e,fi,T 
for charged leptons. Summation over all indices is understood. 



2. Form factors in the MSSM-seesaw 



We present here the analytical results for the form factors in the MSSM-seesaw. 
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where, 
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The couphng factors and self-energies appearing in the neutrahno contributions to the form 
factors are given by, 





— ^Laa^Rab^Rab 


x,X° _ 
''LI — 


x,x° 


— r('*^) n(^) r>{lm)* 

— ^Raa^Lab^Lab 


x,X° _ 
''L2 — 




— ^ Raa-'^ Lab-^ Rab 




^x,x° 


— B^^''^ D^^^ o('m)* 
Laa Rab Lab 




x,x° 


— n^^'^^ n^^) rjCm)* 

~ ^Raa^Rab^Rab 




^x,x° 


nCfc) ri{x) n(lm)* 

— ^ Laa-^ Lab-^ Lab 




^x,x° 


~ ^Laa^Lab^Rab 





^Raa^RPa 
^Laa^LPa 
-^Laa-^RjSa 



9 n /7,.2 „,2 _2 \ D(ifc) Tji^m)* 



X 



(B2) 



The couphng factors and self-energies appearing in the chargino contributions to the form 
factors, K^'^ , L^'^ , and can be obtained from the previous expressions for k^'^° , 
and by making the replacements m^o — > rn^-, mf^ — > rrij)^, B^^^ — > A^^\ D^^'> — > W^^\ 
a ^ i, and b j. 



The form factors -F^ j,, = 1, ..8 can be got from Fj^l., i = 1, ..8 by exchanging the indices 
L ^ R everywhere. 
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